dyn.load("test/kernel.so")

M<-100000
p<-8
n<-200
library(MASS)
X<-as.matrix(Pima.tr[,1:7])
Y<-as.numeric(Pima.tr[,8])
Y[Y==1]<-0
Y[Y==2]<-1
X<-X
v<-apply(X,2,sd)
m<-apply(X,2,mean)
X<-t(apply(t(apply(X,1,"-",m)),1,"/",v))
X<-cbind(1,X)
theta<-matrix(0, nrow=M,ncol=p)

res1<-.C("MHf", as.integer(n),as.integer(p),as.double(Y),as.double(X),as.integer(M))

